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1. Introduction 

In 2005, combined work by the Fermilab Lattice and MILC Collaborations [1] determined the 
value of the D s decay constant fo s to around 10% before it had been determined to that accuracy 
by experiment. When the subsequent experimental determination agreed to within one sigma, 
we claimed that as a successful prediction. As lattice calculations become increasingly accurate, 
of course, at some point we do not expect perfect agreement between the Standard Model and 
experiment. With sufficient precision, the effects of Beyond-the-Standard-Model physics will start 
to show up in low energy measurements. We do not know what that precision will be, so we must 
be cautious in interpreting deviations between theory and experiment. 

Since then, we have increased the precision of our calculations. Our result for fo s remains 
about 10% below the experimental result, and with the increased precision of theory and experi- 
ment, no longer agrees to within one sigma with experiment, as we describe in this paper. Further, 
earlier this year new results on the %, K, D, and D s decay constants appeared from the HPQCD 
Collaboration [2]. They used a new lattice fermion method, Highly Improved Staggered Quarks, or 
"HISQ" fermions, which allowed them to calculate all four decay constants with nearly identical 
methods. They found very good agreement with experiment for the % and K decay constants. Their 
value for the D decay constant was subsequently confirmed by CLEO [3], For fo s , they also found 
a result around 10% below experiment, but with improved precision. Instead of agreement between 
theory and experiment, there is now a greater than three sigma discrepancy. This is the only quan- 
tity in lattice QCD phenomenology with staggered fermions in which such a clear disagreement 
has arisen between theory and experiment, so a puzzle has developed. Precise calculations with 
other lattice methods are of great interest. 

2. Calculations 

We are finishing a reanalysis of the existing data for our calculations of fo, fo s , /b, and fs s 
that is reducing some of our largest uncertainties. We are also preparing for new runs this year with 
four times the statistics. Our calculations are done with improved staggered ("asqtad") light quarks 
[4, 5], and clover/Fermilab [6] &{a) improved heavy quarks. We use the MILC 2 + 1 flavor library 
of unquenched gauge configurations [7], with lattice spacings of around 0.15, 0.12, and 0.09 fm 
(the so-called coarser, coarse, and fine ensembles). The masses of light sea-quarks range between 
0.6m s and 0.lm s . On each of the eleven ensembles, we use from eight to twelve partially quenched 
valence quark masses, ranging from around m s to 0.lm s . 

The decay constants are defined by 
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The combination decay amplitude 
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can be obtained from the correlators 
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The current normalizations are obtained from 

z? = pa r V^^' (2 " 5) 

where Zy® and Z^J are determined nonperturbatively and the remaining (perturbatively calculated) 
short distance corrections in the deviation of from 1 are no more than 0.6%. 
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Figure 1: The leptonic decay amplitudes fV~M for the D and D s mesons, extrapolated to the chiral limit. 
Units are in terms of the heavy-quark potential parameter, r\ . 

Figure 1 shows the extrapolation of the D and D s leptonic decay amplitudes to the physics 
light quark limit. (Units are in terms of the heavy-quark potential parameter r\ .) For ipo (octagons), 
we show only those (fully unquenched) points for which the light valence and sea masses are equal 
to m x , the mass on the abscissa. For (j>D s , we keep both the strange sea mass (mj) and the strange 
valence mass (ms v ) fixed to their simulated values, and plot either as a function of the up/down sea 
mass m.L (crosses), or at ms v (diamonds). The chiral extrapolations make use of all the partially 
quenched data in addition to the points shown. 

The m q dependence is much stronger for the D than the D s , as expected, since in the D s it 
affects only the sea quarks and not the valence quarks. The slope is larger in the continuum limit, 
because taste breaking effects tend to suppress the dependence on the quark mass at finite a. Figure 
2 shows the same thing for the B and B s decay amplitudes, with a qualitatively similar picture. 

Figure 3 shows the extrapolation of the ratio (po / (j>o s to the chiral limit. The slope is strongest 
in the continuum limit (red line and cross), as expected. 

3. Results 

Table 1 shows the uncertainty budgets for the D s , D, B s , and B meson decay amplitudes <j>M, 
and for the ratios Rjj = ^>d/^d, and Rb = <fe/0£,- The three largest uncertainties in our previous 
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Figure 2: The same as for Fig. 1, but for the B and B s mesons. 

results were statistics, heavy quark discretization, and the heavy quark mass. The statistical error 
in the D and D s decay amplitudes has been reduced this year through an improvement in analysis 
method, and without additional data. We are currently incorporating into the chiral and continuum 
extrapolation fits a term for the heavy quark discretization which we expect to substantially reduce 
the uncertainty from this source. The last uncertainty that is large is due to the input heavy quark 
mass, and will be removed with a more careful determination of this quantity. 
We obtain for the decay constants 

f D = 207(11) MeV, (3.1) 

Table 1: Uncertainty budgets in per cent for the leptonic decay amplitudes 0£> s , 0o, 0a s , and and for the 
ratios R D = <fo/<fc> s and R B = <fe/<fe,- 
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Figure 3: The ratio of the D and D s meson leptonic decay amplitudes, extrapolated to the chiral limit. 



In Figure 4, we compare our results for and fy s with the calculations of HPQCD [2] and 
with experiment [3, 8]. For fp, there is very good agreement between experiment, HPQCD, and 
the Fermilab/MILC result. For fn s , there is 

• agreement between HPQCD and Fermilab/MILC, 

• 1.6 a disagreement between Fermilab/MILC and experiment, and 

• 3.5 a disagreement between HPQCD and experiment. 

Many uncertainties cancel in the ratio /d//d s , so we examined this quantity to see if it could 
enhance the significance of the discrepancy between our results and experiment. For now, looking 
at /d//d s doesn't sharpen the picture. In Figure 5, we show our results for /d//d s compared 
with the calculations of HPQCD and with experiment. There is a slight disagreement between 
HPQCD and FNAL/MILC in the ratio, even though and fo s agree within one sigma. Further, 
the experimental uncertainties are independent. They add in quadrature, increasing the size of the 
experimental uncertainty and decreasing the significance of any discrepancy. 




249(11) MeV, 
195(11) MeV, 
243(11) MeV, 



(3.2) 
(3.3) 
(3.4) 



and for the ratios 



f D /f Ds = 0.833(19), 
f B /f Bs = 0.803(28). 



(3.5) 
(3.6) 
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New results for fo and fo s recently appeared from ETMC using twisted-mass fermions [9]. 
They obtained fo = 205 ±10 MeV and fo, = 248 ± 9 MeV, which is in accord with the staggered 
determinations. They present a thorough uncertainty analysis, although we would quibble with 
their use of two rather than three light sea quarks without the inclusion of an uncertainty estimate for 
that approximation. Based on the difference between our unquenched and quenched calculations 
of fo s (249 MeV vs. 213 MeV) [10], we might have guessed a possible uncertainty of 5% from 
leaving out one of the three light sea quarks. (We see charm sea quarks as a different story, since 
m c ~ 1 1 a at our lattice spacings, and the dynamical effects of c quarks are for the most part above 
the cut-off.) 

Three sigma discrepancies between experiment and the Standard Model have occasionally 
appeared and then disappeared before, but the discrepancy in fo, is hard to understand. The un- 
certainty is dominated by experimental statistical error, and three sigma statistical fluctuations are 
very rare. One can double the theory error, and still have a three sigma discrepancy. To explain 
the discrepancy as a theory error, one would have to find a mistake in the theory analysis of fo, 
whose correction would not affect the correct prediction of fo- It is hard to imagine such a mis- 
take. The calculations of fo and fo s are almost identical. The only difference is that fo s should 
be somewhat easier, in that it doesn't require an extrapolation to the physical light quark masses. 
It may be premature to draw ultimate conclusions about the discrepancy, but the result is puzzling 
enough that Kronfeld and Dobrescu have investigated possible new-physics explanations for the 
discrepancy [11]. 

4. Outlook 

We are in the process of reanalyzing our existing data, in which we hope to bring down several 
of our largest uncertainties. New runs are stalling with quadruple the current statistics and at 
smaller lattice spacings, which we expect to help with several of the uncertainties. Comparison of 
fo, in theory and experiment remains a puzzle. This is the only known instance in which lattice 
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Figure 4: Comparison of fo and fo s with the calculations of HPQCD and with experiment. 
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Figure 5: Results for /d//d s compared with the calculations of HPQCD and with experiment. 

QCD with staggered fermions seems to clearly fail to reproduce the Standard Model. This provides 
a good target of opportunity for calculations with other lattice fermion methods. 
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